use replication_data, clear
	
	cap drop occyear
	egen occyear=group(occsoc_num year)
	cap drop occstate
	egen occstate=group(occsoc_num statefip)
	cap drop stateyear
	egen stateyear=group(statefip year)
		
	global covars    ""
	
		capture drop avgageexp
		gen avgageexp=bdshare2000
		capture drop initialneed
		gen initialneed=initial_needbd_plus07
		
		reghdfe dneedbd_plus du  $covars  $weight $restric  , vce(cluster statefip) a(occsoc_num year)
		outreg2 using control_options, replace
				
		reghdfe dneedbd_plus dBGTsupdemratebroad  $covars  $weight $restric  , vce(cluster statefip) a(occsoc_num year)
		outreg2 using control_options, append
		
		*
		reghdfe dneedbd_plus du  $covars  $weight $restric  , vce(cluster statefip) a(occyear)
		outreg2 using control_options, append
				
		reghdfe dneedbd_plus dBGTsupdemratebroad  $covars  $weight $restric  , vce(cluster statefip) a(occyear)
		outreg2 using control_options, append
		*
			
		reghdfe dneedbd_plus du  $covars  $weight $restric  , vce(cluster statefip) a(occyear statefip)
		outreg2 using control_options, append
				
		reghdfe dneedbd_plus dBGTsupdemratebroad  $covars  $weight $restric  , vce(cluster statefip) a(occyear statefip)
		outreg2 using control_options, append
		*
		
			
		reghdfe dneedbd_plus du  $covars  $weight $restric  , vce(cluster statefip) a(occyear occstate)
		outreg2 using control_options, append
				
		reghdfe dneedbd_plus dBGTsupdemratebroad  $covars  $weight $restric  , vce(cluster statefip) a(occyear occstate)
		outreg2 using control_options, append
		
		reghdfe dneedbd_plus du  $covars  $weight $restric  , vce(cluster statefip) a(occyear occstate stateyear)
		outreg2 using control_options, append
				
		reghdfe dneedbd_plus dBGTsupdemratebroad  $covars  $weight $restric  , vce(cluster statefip) a(occyear occstate stateyear)
		outreg2 using control_options, append
			
			
		capture drop avgageexp
		gen avgageexp=avgage2000
		capture drop initialneed
		gen initialneed=initial_need407
		
				reghdfe dneed4 du  $covars  $weight $restric  , vce(cluster statefip) a(occsoc_num year)
		outreg2 using control_options, append
				
		reghdfe dneed4 dBGTsupdemratebroad  $covars  $weight $restric  , vce(cluster statefip) a(occsoc_num year)
		outreg2 using control_options, append
		
		*
		reghdfe dneed4 du  $covars  $weight $restric  , vce(cluster statefip) a(occyear)
		outreg2 using control_options, append
				
		reghdfe dneed4 dBGTsupdemratebroad  $covars  $weight $restric  , vce(cluster statefip) a(occyear)
		outreg2 using control_options, append
		*
			
		reghdfe dneed4 du  $covars  $weight $restric  , vce(cluster statefip) a(occyear statefip)
		outreg2 using control_options, append
				
		reghdfe dneed4 dBGTsupdemratebroad  $covars  $weight $restric  , vce(cluster statefip) a(occyear statefip)
		outreg2 using control_options, append
		*
		
		
			
		reghdfe dneed4 du  $covars  $weight $restric  , vce(cluster statefip) a(occyear occstate)
		outreg2 using control_options, append
				
		reghdfe dneed4 dBGTsupdemratebroad  $covars  $weight $restric  , vce(cluster statefip) a(occyear occstate)
		outreg2 using control_options, append
		
		reghdfe dneed4 du  $covars  $weight $restric  , vce(cluster statefip) a(occyear occstate stateyear)
		outreg2 using control_options, append
				
		reghdfe dneed4 dBGTsupdemratebroad  $covars  $weight $restric  , vce(cluster statefip) a(occyear occstate stateyear)
		outreg2 using control_options, append excel
